Collapse models: analysis of the free particle dynamics 
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We study a model of spontaneous wavefunction collapse for a free quantum particle. We analyze in 
detail the time evolution of the single-Gaussian solution and the double-Gaussian solution, showing 
how the reduction mechanism induces the localization of the wavefunction in space; we also study 
the asymptotic behavior of the general solution. With an appropriate choice for the parameter A 
which sets the strength of the collapse mechanism, we prove that: i) the effects of the reducing 
terms on the dynamics of microscopic systems are negligible, the physical predictions of the model 
being very close to those of standard quantum mechanics; ii) at the macroscopic scale, the model 
reproduces classical mechanics: the wavefunction of the center of mass of a macro-object behaves, 
with high accuracy, like a point moving in space according to Newton's laws. 

PACS numbers: 03.65.Ta, 02.50.Ey, 05.40.-a 



I. INTRODUCTION 

Models of spontaneous wavefunction collapse [J, 0, H, 
1, H i, 0, S I, M EL HI 0, M HI have reached signif- 
icant results in providing a solution to the measurement 
problem of quantum mechanics. This goal is achieved by 
modifying the Schrodingcr equation, adding appropriate 
non-linear stochastic terms 1 : such terms do not modify 
appreciably the standard quantum dynamics of micro- 
scopic systems; at the same time, they rapidly reduce the 
superposition of two or more macroscopically different 
states of a macro-object into one of them; in particular, 
they guarantee that measurements made on microscopic 
systems always have definite outcomes, and with the cor- 
rect quantum probabilities. In this way, collapse models 
describe — within one single dynamical framework - 
both the quantum properties of microscopic systems and 
the classical properties of macroscopic objects, providing 
a unified description of micro- and macro-phenomena. 

In this paper, we investigate the physical properties of 
a collapse model describing the (one-dimensional) evo- 
lution of a free quantum particle subject to spontaneous 
localizations in space; its dynamics is governed by the 
following stochastic differential equation in the Hilbert 
space L 2 (R): 



dip t (x) 



i p 
h 2m 



dl +V\(q- (q) t )dW t 



\{q-(q)t) 2 dt 



(1) 



where q and p are the position and momentum operators, 
respectively, and (q)t = (V'tMV't) denotes the quantum 



average of the operator q; m is the mass of the particle, 
while A is a positive constant which sets the strength 
of the collapse mechanism. The stochastic dynamics is 
governed by a standard Wiener process Wt, defined on 
the probability space (fl, T , P) with the natural filtration 
{Tt : t > 0} defined on it. 

The value of the collapse constant A is given by the 
formula 2 : 



A — — Ao, 
ra 



(2) 



where too is a reference mass which we choose to be equal 
to that of a nucleon, and Ao is a fixed constant which we 
take equal to: 



A 



10 2 m~ 2 sec -1 ; 



(3) 



this value corresponds to the product of the two param- 
eters A GRW and a GRW of the GRW collapse-model 
where A GRW ~ 1CP 16 sec 
a nucleon and l/^/a. 



1 is the localization rate for 
10~ 7 m is the width of the 



Gaussian wavefunction inducing the localizations. 

Eq. ([T]) has already been subject to investigation: in 
ref [lfj a theorem proves the existence and uniqueness 
of strong solutions 3 ; in refs [E fl QU H3] some properties 
of the solutions have been analyzed: in particular, it has 
been shown that Gaussian wavefunctions are solutions of 
Eq. ([T]), that their spread reaches an asymptotic finite 
value (we speak in this case of a "stationary" solution), 
and that the general solution reaches asymptotically a 
stationary Gaussian solution. Finally, in ref p], Eq. ((T|) 
has been first proposed as a universal model of wavefunc- 
tion collapse. 



'Electronic address: bassi@mathematik.uni-muenchen.de 
1 If one aims at reproducing the process of wavefunction collapse 
in measurement-like situation, the terms which have to be added 
to the Schrodinger equation must be non— linear and stochastic, 
since these two are the characteristic features of the quantum 
collapse process. 



2 Differently from [f| , we assume that A is proportional to the mass 
of the particle. 

3 Existence and uniqueness theorems for a wide class of stochastic 
differential equations whose coefficients are bounded operators, 
are studied in ref [l(J . The case of unbounded operators is covered 
in ref |17|1 . See ref [lii . llSil for an introduction to stochastic 
differential equation in infinite dimensional spaces. 
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The aim of our work is to provide a detailed analysis of 
the physical properties of the solutions of Eq. |T]). After 
some mathematical preliminaries (Sec. Hi]) , we will study 
the time evolution of the two most interesting types of 
wavefunctions: the single-Gaussian (Sec. IIII[) and the 
double-Gaussian wavefunctions (Sec. HV[) : in Sec. [V] 
we will discuss the asymptotic behavior of the general 
solution. 

We will next study the effects of the stochastic dynam- 
ics on microscopic systems (Sec. I VI I) and on macroscopic 
objects (Sec. IVII[) ; in the first case, we will see that 
the prediction of the model are very close to those of 
standard quantum mechanics, while in the second case 
we will show that the wavefunction of a macro-object is 
well localized in space and behaves like a point moving 
in space according to the classical laws of motion. We 
end up with some concluding remarks (Sec. IVIII[) . 



II. LINEAR VS NON LINEAR EQUATION 



on (£l,J-, Q), and the Wiener process Wt defined on 
(n,^,P): 

W t = &-2>/A / (q}sds. (7) 
Jo 

The above results imply that one can find the solution 
i\)t of Eq. ([1]), given the initial condition by using the 
following procedure: 

1. Find the solution 4> t of Eq. with the initial 
condition (j) = ip Q . 

2. Normalize the solution: (f> t — > ipt = ^t/W'PtW- 

3. Make the substitution: d£ t -> dW t + 2y/\(q) t . 

The advantage of such an approach is that one can exploit 
the linear character of Eq. ((4]) to analyze the properties 
of the nonlinear Eq. ([1]); as we shall see, the difficult part 
will be computing (q) t from (f> t (step 3): this is where non 
linearity enters in a non-trivial way. 



The easiest way to find solutions of a non-linear equa- 
tion is — when feasible — to linearize it: this is possible 
for Eq. Q and the procedure is well known in the liter- 
ature [1,H, HO, EE EH E3 • Let us consider the following 
linear stochastic differential equation: 



d<t>t(x) 



n 2m 2 



Mx); (4) 



£ t is a standard Wiener process defined on the probability 
space (fi, .F, Q), where Q is a new probability measure, 
whose connection with P will be clear in what follows. 
Contrary to Eq. dJ), the above equation does not pre- 
serve the norm of statevectors, so let us define the nor- 
malized vectors: 



<t>t/\\M if: 
a fixed unit vector if: 



7^0, 

= 0. 



(5) 



By using Ito calculus, it is not difficult to show that ipt 
defined by is a solution of Eq. ([T]) , whenever (f> t solves 
Eq. ([4}. We now briefly explain the relations between 
the two probability measures Q and P, and between the 
two Wiener processes £ t and Wt- 

The key property of Eq. (g]) is that p t = \\4>t\\ 2 is 
a martingale [Tc], E3 satisfying the following stochastic 
differential equation: 



dp t = 2V\(q) t ptd£t 



(6) 



with (q) t — (ipt\q\ipt)- As a consequence of the martin- 
gale property (and assuming, as we shall always do, that 
1 1 (f>o | = 1) pt can be used to generate a new probability 
measure P on (fl, T) [HI]; we choose Q in such a way that 
the new measure P coincides with P. 

Given this, Girsanov's theorem [2l| provides a sim- 
ple relation between the Wiener process £t defined 



III. SINGLE GAUSSIAN SOLUTION 

We start our analysis by taking, as a solution, a single- 
Gaussian wavefunction 4 : 



<j>\ (x) = exp \-a t {x - xtf + ik t x + 



(8) 



where a t and 7 t are supposed to be complex functions of 
time, while x t and k t are taken to be real 5 . By inserting 
© into Eq. one finds the following sets of stochastic 
differential equations 6 : 



da t 

dxt 
dk t 

djf 
dl\ 



A (at) 

m 



dt 

d£t — 2VXxt dt 
-a/X-^I d£t-2V\xtdt 



h — Va r 

™ ktdt+ 2af 
ah 



-772 



Xx 



+Va- 



h —2 

2m ' 



dt + VXxt d£ t - 2\f\x t dt 
dt 



■ x t 



d£ t — 2y\xt dt 



(9) 

(10) 

(11) 
(12) 

(13) 



For a single-Gaussian wavefunction, the two equations 
for j t can be omitted since the real part of j t is absorbed 



* bee ref [H for an analogous discussion within the CSL [j| model 
of wavefunction collapse. 

5 For simplicity, we will assume in the following that the initial 
values of these parameters do not depend on to £ Q. 

6 The superscripts "R" and "I" denote, respectively, the real and 
imaginary parts of the corresponding quantities. 
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into the normalization factor, while the imaginary part 
gives an irrelevant global fase. 

The normalization procedure is trivial, and also the 
Girsanov transformation is easy since, for a Gaussian 
wavefunction like p|). one simply has (q)t — Xt- We then 
have the following set of stochastic differential equations 
for the relevant parameters: 



da t 

dk t = 



A (a t ) 

m 



dt. 



lk t dt+^dW t , 
to 2a" 

-V\^dW t ; 



(14) 

(15) 
(16) 



throughout this section, we will discuss the physical im- 
plications of these equations. 



A. The equation for a t 

Eq. (fll)) for a t can be easily solved [1, [l(| : 
a t = c tanh [b t + k] , 

with: 



(17) 



(1 + 



k 



tanh 



(18) 



After some algebra, one obtains the following analytical 
expressions for the real and the imaginary parts of at : 



R A smh(ujt + ipi) + sin(ajt + ip 2 ) 

' oj cosh(wi + ipi) + cos(ujt + ip2) 

I A sinh(ojt + (f\) — sm(ujt + if2) 

lu cosh(wt + ipi) + cos(u;t + 1P2) ' 

where we have defined the frequency: 



fiA 
mo 



10 



" 5 sec" 1 , 



(19) 
(20) 

(21) 



which does not depend on the mass of the particle. The 
two parameters ipi and (f2 are functions of the initial 
condition: (fx = 2fc R , ip 2 = 2k 1 . 

An important property of af is positivity: 



a" > 



it > V t > 0, 



(22) 



which guarantees that an initially Gaussian wavefunction 
does not diverge at any later time. To prove this, we first 
note that the denominator of (TiT)|) cannot be negative; if 
it is equal to zero then also the numerator is zero: this 
discontinuity can be removed by using expression (fTTl) 
for a t , which — according to the values (flS)) for b and k 
— is analytic for any t. It is then sufficient to show that 
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Schrodinger evolution 




Modified evolution 



hg[t] 



10 



15 



FIG. 1: The picture shows the time evolution (t measured 
in seconds) of the spread in position a q (t) (measured in me- 
ters) of a Gaussian wavefunction, as given by the Schrodinger 
equation and by the stochastic equation. The simulation has 
been made for a nucleon (m = mo). The initial spread has 
been taken equal to 10~ 3 m. 



the numerator remains positive throughout time. Let us 
consider the function f(t) = sxah(u)t + ipi) + sin(wt + 
1P2); we have that /(0) > and f'(t) > for any t, 
which implies that f(t) > for any t, as desired. Note 
that positivity of af matches with the fact that Eq. ((T|) 
preserves the norm of statevectors. 



B. The spread in position and momentum 

The time evolution of the spread in position and mo- 
mentum of the Gaussian wavefunction (|5J| , 




°p® = vV> - (p) 2 = ft > f)2 t K)2 , (23) 

V a t 



is given by the following analytical expressions: 



' cosh(cj£ -(- ipi) + cos(ujt + LP2) 
muj V sinh(wi + <pi) + sin(ujt + (^2) 



a p (t) 



hmuj / cosh(ijt + ipi) — cos(w<: + (P2) 
2 y s'mh(ujt + ipi) + s'm(uit + 1P2) 



(24) 



(25) 



Fig. [T] shows the different time dependence of the spread 
in position, as given by the Schrodinger equation and by 
the stochastic equation: as we see, at the beginning the 
two evolutions almost coincide; as time increases, while 
in the standard quantum case the spread goes to infinity, 
the spread according to our stochastic equation reaches 
the asymptotic value: 



a q (oo) 



H 
muj 




(26) 
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which of course depends on the mass m of the parti- 
cle. This behavior can be understood as follows: the re- 
duction terms (which tend to localize the wavefunction) 
and the standard quantum Hamiltonian (which tends 
to spread out the wavefunction) compete against each 
other until an equilibrium — the stationary solution — 
is reached, which depends on the values of the parameters 
of the model. 

Also the spread in momentum changes in time, reach- 
ing the asymptotic value: 



while the coefficients of covariance matrix 



<Tp (oo) 



hmoj 




Kg m 



sec 



(27) 



It is interesting to compare the two asymptotic values for 
the spread in position and momentum; one has: 



cr 9 (oo)<7 p (oo) = -j= 



(28) 



which corresponds to almost the minimum allowed by 
Heisenberg's uncertainty relations [J, l2Cj : the collapse 
model, then, induces almost the best possible localization 
of the wavefunction — both in position and momentum. 
In accordance with [H, any Gaussian wavefunction hav- 
ing these asymptotic values for a q and a p will be called 
a "stationary solution" 7 of Eq. fl} . 

Note the interesting fact that the evolution of the 
spread in position and momentum is deterministic and 
depends on the noise only indirectly, through the con- 
stant A. 



C. The mean in position and momentum 

The quantities (q) t — x t and (p)t = fikt, corresponding 
to the peak of the Gaussian wavefunction in the position 
and momentum spaces, respectively, satisfy the following 
stochastic differential equations: 



d {q)t = l<p) t<ft + VaJL 



dW t . 



d( P )t 



-Vxn^-dWt. 



(29) 
(30) 



Their average values obey the classical equations 8 for a 
free particle of mass m: 



m±E[(q) t ] = E[(p) t ] 
E[<p)t] = (p>o, 



(31) 
(32) 



7 Strictly speaking, a "stationary" solution is not stationary at all, 
since both the mean in position and the mean in momentum may 
change in time — as it happens in our case; the term "stationary" 
refers only to the spread of the wavefunction. 

8 These equations can be considered as the stochastic extension of 
Ehrenfest's theorem. 



C(t) 



E 



(q)t-E[(q) t ] 

(p)t-n( P )t] 

Cg*(t) C qp (t) 

C pq (t) C p 2(t) 



(ff>t-E[<9>t] 

<p)t-n<p) t ] 



evolve as follows 

d_ 
dt 
d_ 
~dt 

d 



C qp (t) 



dt 



C p 2 (t) 



rn 
1 

m 

Xh 2 



C p 2 (t) 



4 (at) 2 ' 
Xh a\ 

TJaW' 



(33) 
(34) 

(35) 



Particularly interesting is the third equation, which im- 
plies that the wavefunction picks larger and larger com- 
ponents in momentum, as time increases; as a conse- 
quence, the energy of the system increases in time, as 
it can be seen by writing down the stochastic differential 
for (H) t = (p 2 ^' 



d(H) t 



/t/2m: 

_ Xtf 
2m 



dt - a/A - 



(p)tdW t , (36) 



from which it follows that: 

dt 2m 2mo 



10 



-43 



J/sec. (37) 



This energy non-conservation is a typical feature of 
space-collapse models, but with our choice for the pa- 
rameter A, the increase is so weak that it cannot be de- 
tected with present-day technology (see ref. [l|). 



IV. DOUBLE GAUSSIAN SOLUTION 

We now study the time evolution of the superposition 
of two Gaussian wavefunctions; such an analysis is inter- 
esting since it allows to understand in a quite simple and 
clear way how the reduction mechanism works, i.e. how 
the superposition of two different position states is re- 
duced into one of them. To this purpose, let us consider 
the following wavefunction: 

4>°{x) = 4>xt(x) + (f>2t(x) = 

= exp [-a u (x - xitf + ik lt x + j u ] 

+ exp [-a 2t (x - X2t) 2 + ikitx + 72*] ; (38) 

we follow the strategy outlined in Sec. II, by first finding 
the solution of the linear equation. 

Because of linearity, 4>f (x) is automatically a solution 
of Eq. (|4]), provided that its parameters satisfy Eqs. ([9]) 
to (Tj~3")) . The difficult part of the analysis is related to 
the change of measure: the reason is that in the double- 
Gaussian case the quantum average (q)t is not simply 
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equal to xu or X2t, as it is for a single-Gaussian wave- 
function, but is a nontrivial function 9 of all the param- 
eters denning <f)f (x) ; in spite of this difficulty, the most 
interesting properties of the dynamical evolution of 4>f (x) 
can be analyzed in a rigorous way. 

We first observe that the two equations for an and an 
are deterministic and thus insensitive of the change of 
measure; accordingly, the spread (both in position and 
in momentum) of the two Gaussian functions defining 
4>t{x) evolve independently of each other, and maintain 
all the properties discussed in the previous section. For 
simplicity, we assume that ait — &2t at t — so that these 
two parameters will remain equal at any subsequent time. 



A. The asymptotic behavior 

Let us consider the differences Xt — X2t — xu and 
Kt — ^2t ~ kit between the peaks of the two Gaussian 
functions in the position and in the momentum spaces; 
they satisfy the following set of equations: 



' X t ' 









-A x {t) h/m 
-A 2 {t) 





' x t ' 







dt. 



with: 



Mt) 



= 2X 



cosh(u/t + ipi) + cos(ujt + tp 2 ) 
sinh(a;£ + ipi) + sin(u;i + ip%) 
sinh(o;< + fi) — sin(u;t + ip 2 ) 
sinh(wt + t/?i ) + sin(o>i + ^2) 



(39) 

(40) 
(41) 



We see — this is the reason why we have taken into ac- 
count the differences X t and K t — that the above system 
of equations is deterministic, so it does not depend on the 
change of measure. 

The coefficients of the 2x2 matrix A(t) defining the 
linear system ([3TR) are analytic in the variable t, and the 
Liapunov's type numbers [24| of the system are the same 
as those of the linear system obtained by replacing A(t) 
with A(po), where: 



A(oo) 



lim A(t) 

t — >+oo 



— lu h/m 
-2A 



The eigenvalues of the matrix A(oo) are: 

1 

2 



Mi, 2 



(l±i) w, 



(42) 



(43) 



from which it follows that the linear system (|39[) has only 
one Lyapunov's type number: —u>/2. This implies that, 
for any non trivial (vector) solution Z(t) of (|3"9"|) . one has: 



lim 

t — >+oo 



ln|Z(t)| 



u> 
2 



lim Z(t) 



0. (44) 



9 This is the reason why the dynamics of these parameters changes 
in a radical way, with respect to the single— Gaussian case. 



We arrive at the following result: asymptotically, the dif- 
ference X2t — xu between the peaks of the two Gaussian 
wavefunctions in the position space goes to zero; also the 
difference k2t — kit between the peaks of the two Gaus- 
sian wavefunctions in the momentum space vanishes. In 
other words, the two Gaussian wavefunctions converge 
toward each other, and asymptotically they become one 
single-Gaussian wavefunction which, from the analysis of 
the previous subsection, is a "stationary" solution of the 
stochastic equation |T|). 

Anyway, this behavior in general cannot be responsible 
for the collapse of a macroscopic superposition; as a mat- 
ter of fact, let us consider the following situation which, 
for later convenience, we call "situation A" : 

i) The wavefunction is in a superposition of the form (|38|) 
and, at time t = (consequently, also for any later t), 
an and a 2 t are equal to their asymptotic value; 

ii) The distance K t = k%t ~ kit is zero at time t = 0. 
Under these assumptions, the linear system (|39p can be 
easily solved and one gets: 



Xi 



X e 



-ut/2 



tot 



(45) 



We see that the time evolution of X t is independent of 
the mass of the particle: this result implies that, when 
the spread of the two Gaussian wavefunctions is equal 
to its asymptotic value 10 , their distance decreases with 
a rate oj/2 ~ 10~ 5 sec -1 , which is too slow to justify a 
possible collapse, in particular at the macro-level. 



B. The collapse 

Now we show that the collapse of the wavefunction oc- 
curs because, during the evolution, one of the two Gaus- 
sian wavefunctions is suppressed with respect to the other 
one 11 ; the quantity which measure this damping is the 
difference Tf — 7^ — 7f t , which satisfies the stochastic 
differential equation: 



dTf = -XX t [x u + x 2t - 2(q) t ] dt + VXXtdWt; (46) 

if Tf — > +00, then <pi is suppressed with respect to 4>2 and 
the superposition practically reduces to <fi 2 ; the opposite 
happens if Tf — * —00. 

To be more precise, we introduce a positive constant 
b which we assume to be conveniently large (let us say, 
b = 10) and we say that the superposition is suppressed 
when |r R | > 6; moreover, we say that: 



6 D is reduced to 
b° is reduced to 



>i when: 
'~>2 when: 



pa < 

r R > 



10 We will see in Sec. IVIII that at the macroscopic level the spread 
of a Gaussian wavefunction converges very rapidly towards its 
asymptotic value. 

11 Accordingly, the collapse mechanism is precisely the same as the 
one of the original GRW model [l| . 
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We now study the time evolution of Tf. 

By writing (q) t in terms of the coefficients defining 



2a? 



i-) R. n R, R | R 



(47) 



it is not difficult to prove that Eq. (|4"6")> becomes: 

dTf = XX 2 tanhrf dt + g t (T?)dt + VxX t dW t , (48) 
where we have defined the following quantities: 



St 



and: 

fh = 

Y t = 

t = 
with a? 



h t [(xit + x 2 t) cos f t 
2XX t h f e r * ■ 



6 t +Yt sin t 



(49) 
(50) 



sin ft 



1 



„2r? 



+ X t cos ( 



[f + e 2r * + 2ht cos 9 t e r ?] [l + e 2r *] 



exp 



2a\X t + K t 
2a? ' 



'2/ • 



(51) 
(52) 
(53) 



Eq. (|48|) cannot be solved exactly, due to the presence 
of the term g t (Tf) which is a non-simple function of T?; 
to circumvent this problem, we proceed as follows. We 
study the following stochastic differential equation: 



of? = XX? tanh f? eft + VxX t dW t , 



(54) 



which corresponds to Eq. (|48|) without the term g t (T?), 
and at the end of the subsection we estimate the error 
made by ignoring such a term. 

In studying Eq. (|5"4"|) , it is convenient to introduce the 
following time-change: 



t 



s t = X XUu. 



(55) 



X 2 is a continuous, differentiable and non-negative func- 
tion, which we can assume not to be identically zero 12 in 
any sub-interval of ]R + ; as a consequence, St is a mono- 
tone increasing — thus invertible — function of t and Eq. 
([55)1 defines a good time change. 

Under this time substitution, Eq. (|5"4"|) becomes: 



df* = tanhf^ds + dW s , 



(56) 



12 If there exists an interval / of K + such that Xt = Vt £ /, then 
Eq. I|39j l implies that Xt remains equal to for any subsequent 
time, i.e. the two Gaussian wavefunctions coincide. 



where 

W s = \/A f X t dWt (57) 
Jo 

is a Wiener process 13 with respect to the filtered space 
{ft, J 7 , {T s : s > 0},P). Note that, since according to 
Eq. (14"4"1) , X t in general decays exponentially in time, 
Soo < oo, and Eq. (|56|) is physically meaningful only 
within the interval 14 [0, Sqo). 

Eq. (|56p can be analyzed in great detail [25[; in par- 
ticular, the following properties can be proven to hold 15 : 
1. Let us define the collapse time Sb = inf{s : \T?\ > b}: 
this time is finite with probability 1 and its average value 
is equal to 16 : 



E[sy 



b tanh b — bo tanh h 



(58) 



If bo <C b (which occurs when both terms of the superpo- 
sition give a non-vanishing contribution), and since we 
have assumed b ~ 10, then E[5b] ~ b. 
2. The variance Y[S b ] = E[S 2 ] - E 2 [Sb] is given by 17 : 



F(b) - F(bo), 



with: 



F (x) — x 2 tanh 2 x + x tanh x — x 2 



(59) 



(60) 



Note that since F(x) is an even, positive function, in- 
creasing for positive values of x, it follows that Y[Sb] > 
whenever |6 | < b, as it assumed to be. 
3. The collapse probability P^ 2 that </> D is reduced to <f>2, 
i.e. that hits point b before point —6 is given by 18 : 



1 tanh b + tanh bo 

2 tanh b 



(61) 



according to our choice for b, tanh& ~ I, and conse- 
quently 



-[l+tanh&n] 



'10 



20 



(62) 



which (neglecting the overlapping between the two Gaus- 
sian wavefunctions) corresponds to the standard quantum 
prescription for the probability that cjf collapses to 02- 



13 See pages 111-113 of [2Sl . 

14 Thus, strictly speaking, Eq. J5TH defines a Wiener process only 
for s < s^; we then extend in a standard way the process to the 
interval (soo , oo). 

15 Throughout the analysis, we will assume that Tq = bo £ (—6, 6). 

2l. 
251. 



16 See theorem 2, page 108 of 

17 See theorem 3, page 109 of 
18 



See theorem 4, page 110 of 25 
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4. The delocalization probability, i.e. the probability P^ 1 

that Tf goes below b — r\ before time Sqc, after having 
reached b (77 is a positive quantity smaller than 6), is 19 : 



^>2 



< 1 - P{inff* >b-rr} = 

, , . tanh 77 
= 1 - (1 + tanhb) ' 



1 + tanh i] 



(63) 



In the above definition, we have required a delocalization 
to occur before time Sqo because — since Soc corresponds 
to t = oo — a delocalization at a time > does not 
correspond to a real physical delocalization. 

If, for example, we take r\ = 3, we have that P? el < 
0.002. 

This concludes our analysis of the statistical behavior 
of r^. We end this section by discussing how one can 
estimate the error made in neglecting the term g t (Tf) in 
Eq. (|48]) . i.e. in studying in place of T™. In Appendix 
A the following estimate for gt(Tf) is given: 



\9t(Tf)\ < A 



\X t \ + \Y t \f 



< 



e af(\X t \ + \Y t \)y4 

xxf 



(64) 
(65) 



are the minimum values af and X t 



where and X 
take during the time interval one is considering. For 
example, if we take a 1-g object and we assume that 

a. 



and X m — 1 cm, then we have: 



1 



(66) 



which is very close to zero. 

By using inequality (|6"5|) and lemma 4, pag. 120 of [25| . 
one can easily show that T^~~ < < r*+, where 
are solutions of the following two stochastic differential 
equations: 



dF^ = [tanhT^ic] ds 



dW K 



(67) 



with an obvious meaning of the signs (as before, we have 
moved from the variable t to the variable s). 

Eqs. (f6"T)) can be analyzed along the same lines fol- 
lowed in studying Eq. (f56|) . getting basically the same 
results, due to the very small value c takes in most rel- 
evant physical situations. This kind of analysis will be 
done in detail in a future paper: there, we will study 
the most important case where a macroscopic superpo- 
sition can be created, i.e. a measurement -like situation 
in which a macroscopic object acting like a measuring 
apparatus interacts with a microscopic system being ini- 
tially in a superposition of two eigenstates of the operator 
which is measured; we will show that, throughout the in- 
teraction, the wavefunction of the apparatus is — with 



extremely high probability — always localized in space, 
and that the measurement has a definite outcome with 
the correct quantum probabilities. 



V. GENERAL SOLUTION: THE ASYMPTOTIC 
BEHAVIOR 

In this section we analyze the asymptotic behavior of 
the general solution of Eq. ([T]), showing that — as time 
increases — any wavefunction collapses towards a sta- 
tionary Gaussian solution. 



A. The collapse 

We have seen in the previous sections that both the 
single-Gaussian and double-Gaussian solutions asymp- 
totically converge towards a stationary solution having 
the form: 



with: 



2aJL 



exp 



lim a t 

t — > + oo 



(x-x t ) 2 + ik t x] , (68) 



(69) 



it becomes then natural to ask whether such a kind of 
wavefunction is the asymptotic limit of any initial wave- 
function. The answer is positive [2(| as we shall now see 
by following the same strategy used in refs. [ll|, Q H3| 
to prove convergence of solutions. 

Since a wavefunction of the form (|68p is an eigenstate 
of the operator: 



A 



1 



■li- 



the proof consists in showing that the variance 20 : 
AA t ^ & t \[Ai -(Ai)}[A-(A)U t ) 



(70) 



(71) 



converges to for t 
AA t holds: 



-oo. The following expression for 



AA t 
with: 



Aq t 



-Apt 



-^{q,p) 



h 



(72) 



S(g,p) 



[<M«-<«>][p-<p>]W*> 



(73) 



See formula 6, page 117 of [25 



20 In fact, it is easy to prove that AAt = if and only if \ipt) is an 
eigenstate of the operator A, from which it follows that AAt — > 
if and only if |i/>t) converges towards an eigenstate of A. In Eq. 
frO , (A) = (i()t\Aipt), and similarly for (At). 
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After a rather long calculation, one finds that: 

d_ 

di 



'' ^[AA t ] = -ELAA t + 2A (Aq t - a 2 q (oo)) 2 



-2A Aq 



^(q,p) 



< 0. 



(74) 



Since, by definition, E [AAj] is a non-negative quantity, 
Eq. (|74|) is consistent if and only if the right-hand-side 
asymptotically vanishes. This, in particular, implies that 



AA t 



0. 



(75) 



t — > + oo 

except for a subset of f2 of measure zero. 

B. The collapse probability 

We now analyze the probability that the wavefunction, 
as a result of the collapse process, lies within a given re- 
gion of space 21 . To this purpose, let us consider the fol- 



lowing probability measure, defined on the Borel sigma- 
algebra B(R): 



IH(A) = E P [\\Pa1h 



(76) 



where Pa is the projection operator associated to the 
Borel subset A of M. Such a measure is identified by the 
density p t (x) = E v [\ip t {x)\ 2 ]: 



Pt(A) = / p t (x)dx, 

J A 



(77) 



and it can be easily shown that the density pt(x) cor- 
responds to the diagonal element (x\pt\x) of the sta- 
tistical operator 22 ^ = Ep [|-0t)(^t|], which satisfies the 
Lindblad-type equation: 



dt Pt 



■tkh (78) 



The solution of the above equation, expressed in terms of 
the solution (qi\ps{t)\Q2) of the pure Schrodinger equa- 
tion (A = 0), is 23 : 



i r+oo r+oc 

<91|P(*)I«2> = TTT / dk / dye-^ k yF[k 1 q 1 -q 2 ,t](q 1 +y\p s (t)\ q2 +y), 



(79) 



with: 



F[k, x,t] — exp 



A / 2 k k 2 , 

2 V m 3m 2 



(80) 



r 



The Hermitian symmetry of (qi\p(t)\q2) follows from the 
fact that F[k,x,t] = F[—k,—x,t}; for A = we have 
F[k,X,t] = 1 so that {qi\p(t)\q2) = (qi\ps(t)\q2) as it 
must be. 

From Eq. ([75]) it follows that: 

i r+oo r+oc 

P^ X ) = ^ dk dye-^ k yF[k,0,t]pf(x + y) 

(81) 

where pf (x) = (x\ps(t)\x) is the standard quantum prob- 
ability density of finding the particle located i in a po- 
sition measurement. One can easily perform the integra- 
tion over the k variable and he gets: 



«-t-oo 



dye-^y 2 pf(x + y), 



For a macroscopic object (let us say m > 1 g) and for 
very long times (e.g. t £ [0, 10 8 sec ~ 3 years], a time 
interval which of course is much longer than the time 
during which an object can be kept isolated so that the 
free particle approximation holds true) at is a very large 
number, so large that the exponential function in Eq. (1821) 
is significantly more narrow than pf(x), for most typical 
wavefunctions (in Sec. lVII Al we will see that the asymp- 
totic spread of the wavefunction of a 1-g object is about 
10 -13 m). Accordingly, this exponential function acts 
like a Dirac-delta and one has: pt(x) ~ pf (x), with very 
high accuracy. This in turn implies that: 



MA) c \\P A ^ 



Sch||2 



(84) 



(82) 



with: 



a i 



3toq to ^ 43 / m\ /sec y 
2h 2 A i 3 ~ I kg J \ t ) 



(83) 



21 Of course, wavefunctions in general do not have a compact sup- 
port; accordingly, saying that a wavefunction lies within a given 
region of space amounts to saying that it is almost entirely con- 
tained within the region, except for small "tails" spreading out 
of that region. 
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in other words, the probability measure /it (A) is very 
close to the quantum probability of finding the particle 
lying in A in a position measurement. 

We still have to discuss the physical meaning of the 
probability measure /it (A) defined by Eq. (|76|) : such a 
discussion is relevant only at the macroscopic level, since 
we need only macro-objects to be well localized in space 
(and with the correct quantum probabilities). 

As already anticipated the spread of a wavefunction of 
a macro-object having the mass of e.g. 1 g reaches in 
a very short time a value which is close to the asymp- 
totic spread <7 g (oo) ~ 10~ 13 m; then, if we take for A an 
interval [a, b] whose length is much greater than cr g (oo) 
- for example, we can take b — a = 10~ 7 m, which is 
sufficiently small for all practical purposes — only those 
wavefunctions whose mean lies around A give a non van- 
ishing contribution to /it (A). As a consequence, with the 
above choices for A (and, of course, waiting a time suffi- 
ciently long in order for the reduction to have occurred) , 
the measure /it (A) represents a good probability measure 
that the wavefunction collapses within A. 



Cause of decoherence 


10" 3 cm 10" 6 cm 
dust particle large molecule 


Air molecules 
Laboratory vacuum 
Sunlight on earth 
300K photons 
Cosmic background rad. 


10 36 10 30 
10 23 1 17 
10 21 10 13 
10 19 10 6 
10 6 KT 12 


COLLAPSE 


10 7 1(T 2 



TABLE I: A in cm" sec" for decoherence arising from dif- 
ferent kinds of scattering processes (taken from Joos and Zeh 
[2^]L In the last line: A for the collapse model given by Eq. 
©■ 



2. The spread of the wavefunction reaches an asymptotic 
value which depends on the mass of the particle: 



a q (oo) 
(Tq(oo) 



I m for an electron, 
1 cm for a nucleon. 



VI. EFFECT OF THE REDUCING TERMS ON 
THE MICROSCOPIC DYNAMICS 

In the previous sections we have studied some analyti- 
cal properties of the solutions of the stochastic differential 
equation ([T]); we now focus our attention on the dynamics 
for a microscopic particle. 

The time evolution of a (free) quantum particle has 
three characteristic features: 

1. The wavefunction is subject to a localization pro- 
cess, which, at the micro-level, is extremely slow, almost 
negligible. For example, with reference to the situation 
A we have defined at the end of Sec. IIVA1 Eq. fl5j) 
shows that the distance X t between the centers of the two 
Gaussian wavefunctions remains practically unaltered for 
about 10 sec; under the approximation X t ~ Xq, Eqs. 
([55)) and ([55]) imply that the time necessary for one of 
the two Gaussians to be suppressed is: 

/ m \ 2 

E [Xy — 10 6 I — I sec for an electron, 
\ x oJ 

E [Xy ~ 10 3 (j^J sec for a nucleon, 

which are very long times 24 , compared with the charac- 
teristic times of a quantum experiment. 



The definition pt = Ep [\ipt) (ipt |] is rather formal; see [161I17H for 
a rigorous definition. 

Ref. pj|, Appendix C, shows how the solution can be obtained. 
In the first case, indeed, the reduction time is longer than the 
time during which the approximation Xt ~ Xq is valid — see 
Eq. {45}. 



3. The wavefunction undergoes a random motion both 
in position and momentum, under the influence of the 
stochastic process. These fluctuations can be quite rele- 
vant at the microscopic level; for example, for a station- 
ary solution, one has from Eq. (|33|) : 

-j-V[(c/)t] > 1 m 2 /sec for an electron, 
^V[(g)t] > 10" 3 m 2 /sec for a nucleon. 

This is the physical picture of a microscopic particle 
as it emerges from the collapse model. In order for the 
model to be physically consistent, it must reproduce at 
the microscopic level the predictions of standard quan- 
tum mechanics, an issue which we are now going to dis- 
cuss. 

Within the collapse model, measurable quantities are 
given [? ] by averages of the form Ep[(0) t ], where O is 
(in principle) any self adjoint operator. It is not difficult 
to prove that: 

E P [(0) t ] = Tr [Op t ], (85) 

where the statistical operator p t = Ep [|V't)(V't|] satis- 
fies the Lindblad-type Eq. ([78]) . This is a typical master- 
equation used in decoherence theory to describe the inter- 
action between a quantum particle and the surrounding 
environment (26j ; consequently, as far as experimental re- 
sults are concerned, the predictions of our model are sim- 
ilar to those of decoherence models 25 . It becomes then 



We recall the conceptual difference between collapse models and 
decoherence models. Within collapse models, one modifies quan- 
tum mechanics (by adding appropriate nonlinear and stochastic 
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natural to compare the strength of the collapse mecha- 
nism (measured by the parameter A) with that of deco- 
herence. 

Such a comparison is given in Table HI when the system 
under study is a very small particle like an electron, or 
an almost macroscopic object like a dust particle. We see 
that, for most sources of decoherence, the experimentally 
testable effects of the collapse mechanism are weaker than 
those produced by the interaction with the surrounding 
environment. This implies that, in order to test such 
effects, one has to isolate a quantum system — for a suf- 
ficiently long time — from almost all conceivable sources 
of decoherence, which is quite difficult: the experimen- 
tally testable differences between our collapse model and 
standard quantum mechanics are so small that they can- 
not be detected unless very sophisticated experiments are 
performed [HHU. 



VII. MULTI-PARTICLE SYSTEMS: EFFECT OF 
THE REDUCING TERMS ON THE 
MACROSCOPIC DYNAMICS 

The generalization of Eq. (Q]) to a system of N inter- 
acting and distinguishable particles is straightforward: 



d^ t {{x}) = 



N 



-H T dt + J2 (q n - (q n )t) dW r t 



1 N 



Mi*}), (86) 



where H T is the standard quantum Hamiltonian of the 
composite system, the operators q n (n = 1, . . . N) are the 
position operators of the particles of the system, and 
(n = 1, . . . N) are N independent standard Wiener pro- 
cesses; the symbol {x} denotes the N spatial coordinates 
xi, . . .x N . 

For the purposes of our analysis, it is convenient 
to switch to the center-of-mass (R) and relative 
(xi, X2, • ■ ■ xn) coordinates: 



R 



- N N 

-tt ^m„x„ M = y^m» 



71=1 



71=1 



(87) 



R ~t~ x n , 



terms) so that macro-objects are always localized in space. De- 
coherence models, on the other hand, are quantum mechanical 
models applied to the study of open quantum system; since they 
assume the validity of the Schrodinger equation, they cannot in- 
duce the collapse of the wavefunction of macroscopic systems (as 
it has been shown, e.g. in [2?|l even if one of the effects of the 
interaction with the environment is to hide — not to eliminate 
— macroscopic superpositions in measurement-like situations. 



let Q be the position operator for the center of mass and 
q n (n = 1 . . . N) the position operators associated to the 
relative coordinates. 

It is not difficult to show that — under the assumption 
H T = H cu + H Tel — the dynamics for the center of mass 
and that for the relative motion decouple; in other words, 
ip t ({x}) = ipt™( R ) ® iPfdx}) solves Eq. ([86]) whenever 
ipt™(R) and ^ o1 ({x}) satisfy the following equations: 



N 



-- H le} dt + J2 VKtin - {q n )t)dW t 



1 N 



dt 



iPf({x}),(88) 



with: 



- T H C m dt - 
n 



(Q - (Q)t) 2 dt 



N 



— \i — — Ao. 

^ mo 

n— 1 



(89) 



(90) 



The first of the above equations describes the internal 
motion of the composite system, and will not be analyzed 
in this paper; in the remainder of the section, we will 
focus our attention on the second equation. 

Eq. (|89|) shows that the reducing terms associated to 
the center of mass of a composite system are equal to 
those associated to a particle having mass equal to the to- 
tal mass M of the composite system; in particular, when 
the system is isolated — i.e., H CM = P 2 /2M, where P is 
the total momentum — the center of mass behaves like a 
free particle, whose dynamics has been already analyzed 
in Sects. IIII1 ITVl andlVl In the next subsections we will 
see that, because of the large mass of a macro-object, 
the dynamics of its center of mass is radically different 
from that of a microscopic particle. 



A. The amplification mechanism 

The first important feature of collapse models is what 
has been called the "amplification mechanism" [l], [3] : the 
reduction rates of the constituents of a macro-object sum 
up, so that the reduction rate associated to its center 
of mass is much greater than the reduction rates of the 
single constituents. 

This situation is exemplified in Fig. [21 which shows the 
time evolution of the spread <J q {i) of a Gaussian wave- 
function. We note that, however large the initial wave- 
function is, after less than 10 -2 sec — which corresponds 
to the perception time of a human being — its spread 
goes below 10~ 5 cm, which is the threshold chosen in 
the original GRW model [l[ , below which a wavefunction 
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FIG. 2: Time evolution of the spread <J q (t) of a Gaussian 
wavefunction of the center of mass of a system containing 
N = 10 24 nucleons; different initial conditions have been con- 
sidered. Time is measured in seconds, while the spread is 
measured in meters. 



can be considered sufficiently well localized to describe 
the classical behavior of a macroscopic system. 
More generally, Eq. If71]) implies that: 



E [A A t ] 



> 2A [Aq - a 2 q (^)Y 



(91) 



and, as long as the spread of the wavefunction is signif- 
icantly greater than its asymptotic value, i.e. the wave- 
function is not already sufficiently well localized in space, 
we have: 



— E [AA ( ] 



'»=■' ' i) (£) » 2 /»c. OB) 



which is a very high reduction rate, for a macroscopic 
object. Note that, as previously stated, the velocity in- 
creases, for increasing values of the mass of the particle. 

Asymptotically, the wavefunction of a macro-object 
has an extremely small spread; for example: 



(Tg(oo) 
a q (oo) 



1CT 13 m 



10 



-27 



for an 1-g object, 
for the Earth. 



Thus, according to our collapse model, macro-particles 
behave like point-like particles. 



B. Damping of fluctuations 

We have seen that the mean (both in position and in 
momentum) of the wavefunction undergoes a diffusion 
process arising from the stochastic dynamics: such a dif- 
fusion is quite important at the microscopic level, and 
it is responsible for the agreement of the physical predic- 
tions of the model with those given by standard quantum 
mechanics. We now analyze the magnitude of the fluctu- 
ations at the macroscopic level. 



Contrary to the behavior of the reduction mechanism, 
which is amplified when moving from the micro- to the 
macro-level, the fluctuations associated to the motion 
of microscopic particles interfere destructively with each 
other, in such a way that the diffusion process associated 
to the center of mass of an iV-particle system is much 
weaker than those of the single components. We now 
give some estimates. 

Let us suppose that the center-of-mass wavefunction 
has reached a stationary solution; under this assumption, 
one has from Eqs. ([33]), (J34J) and ([35]): 



V[(p) t ] - xnh. 



{tot) 3 {cot) 2 



3! 



2! 



tot 



(93) 
(94) 



Since, for example, oj/8X ~ 10~ 27 m 2 for a 1-g object, 
and uj/8X ~ 10~ 54 m 2 for the Earth, we see that for a 
macro-object the numerical values of the parameters are 
such that for very long times (in many cases much longer 
that the age of the universe) the fluctuations are so small 
that, for all practical purposes, they can be neglected; 
this is how classical determinism is recovered within our 
stochastic model. 

The above results imply that the actual values of (q)t 
and (p)t are practically equivalent to their stochastic av- 
erages; since these stochastic averages obey the classical 
laws of motion ([3T]) and ([32]) . we find out that (q)t and 
(p)t practically evolve according to the classical laws of 
motion, for most realizations of the stochastic process. 

The conclusion is the following: in the macroscopic 
regime, the wavefunction of a macroscopic system be- 
haves, for all practical purposes, like a point-like parti- 
cle moving deterministically according to Newton laws of 
motion. 



VIII. CONCLUSIONS 

From the analysis of the previous sections we have seen 
that, in general, the evolution of the wavefunction as pre- 
dicted by the collapse model is significantly different from 
that predicted by standard quantum mechanics, both at 
the micro- and at the macro-level. For example, at the 
microscopic level the random fluctuations can be very 
large, while in the standard case there are no fluctua- 
tions; at the macroscopic level, wavefunctions rapidly lo- 
calize in space, while in the standard quantum case they 
keep spreading. 

Anyway, as far as physical predictions are concerned, 
our model is almost equivalent to standard quantum me- 
chanics, the differences being so small that they can 
hardly be detected with present-day technology. More- 
over, at the macroscopic level the localization mechanism 
becomes very rapid and the fluctuations almost disap- 
pear: the wavefunction of the center of mass of a macro- 
scopic object behaves like a point-like particle moving 
according to Newton's laws. 
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To conclude, the stochastic model reproduces, with ex- 
cellent accuracy, both quantum mechanics at the micro- 
scopic level and classical mechanics at the macroscopic 
one, and describes also the transition from the quantum 
to the classical domain. 
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is bounded between 1 and —1, we get: 



\9t(T?)\ < A 



hi 



< X 



1 - ht 

ht 
1-ht 



[X 2 + \X t Y t \] (A4) 



l* t | + |y t | 
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APPENDIX A: DERIVATION OF 
INEQUALITY ([64]) 

From Eq. (|49[) . one has: 

Is* 0T) I < 



(Al) 



2X\X t \ht 



\X t cos Of, 



lit sin t 



l + e 2r ? + 2h t cos0 t e r t 



e * 



l + e 2r * + 2h t cos6» t e r ? 



1-e 



21, 



1 + e 



2r? 



By using the fact that sin# t < 1 and cos Ot < 1, that the 
function 



x>0, -1<c<1 (A2) 



1 + x 2 + 2cx 

is bounded between and 1/(2 + 2c), while the function 

1 - x 



l+x 



x>0 



(A3) 



from which, dividing both the numerator and the 
denumerator on the right-hand-side by h t and using 
(x + y) 2 < 2(x 2 + y 2 ), inequality (f6"4")l follows. 

NOTE (FIRST PUBLISHED IN: A. BASSI, D.G.M. 
SALVETTI, J. PHYS. A: 40, 9859 (2007)) 



Eqs. (|12|) and (| 13[) contain a computational mistake. 
The correct equations are: 



\xi 



h . A 



■in 



dt 



V\x t d£ t — 2y/Xx t dt 



h „ h Xal 
m 1 2m 1 4(of ) 2 
a 1 r 

+V\^-x t d£ t -2V\x t dt 



dt 



(A5) 



(A6) 



which differ from Eqs. (|12|l and (|13ll . in the first case 
for the extra factor A/4af and in the second case for 
the factor Xa I t /A(af) 2 . We correct in this way a mistake, 
which however does not affect the subsequent analysis. 
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